Paramagnetic ground state with field-induced partial order in Nd3Ga5SiOi4 probed by 

low-temperature heat transport 



(N 

O 
(N 

>>: 



u 



o 



rn 

in 
o 



Q. J. Li,i Z. Y. Zhao,! H. D. Zhou,^ W. P. Ke,i X. M. Wang i 
C. Fan,i X. G. Lin,^ L. M. Chen,^ X. Zhao,4 and X. F. Sun^'Q 

' Hefei National Laboratory for Physical Sciences at Microscale, 
University of Science and Technology of China, Hefet, Anhm 230026, People's Republic of China 
^National High Magnetic Field Laboratory, Florida State University, Tallahassee, Florida 32306-4005, USA 
^Department of Physics, University of Science and Technology of China, 
Hefei, Anhui 230026, People's Republic of China 
^School of Physical Sciences, University of Science and Technology of China, 
Hefei, Anhui 230026, People's Republic of China 
(Dated: March 4, 2013) 

We study the low-temperature heat transport of Nd3Ga5SiOi4, which is a spin-hquid candidate, 
to probe the nature of ground state and the effect of magnetic field on the magnetic properties. The 
thermal conductivity («:) shows a purely phononic transport in zero field. The external magnetic 
field along the c axis induces a dip- like behavior of k{H), which can be attributed to a simple 
paramagnetic scattering on phonons. ffowever, the magnetic field along the ab plane induces another 
step-like decrease of k. This kind of n{H) behavior is discussed to be related to a field-induced partial 
order, which yields low-energy magnetic excitations that significantly scatter phonons. These results 
point to a paramagnetic ground state that partial magnetic order can be induced by magnetic field 
along the ab plane, which is also signified by the low-T specific heat data. 

PACS numbers: 66.70.-f, 75.47.-m, 75.50.-y 



I. INTRODUCTION 

Geometrically frustrated magnetic materials with com- 
peting spin interactions give rise to macroscopic degen- 
eracy of ground state. They have displayed abundant 
exotic magnetic phenomena and introduced important 
concepts to condensed matter physics The recently 
discovered rare-earth Langasites, i?3Ga5SiOi4 (i? — Nd, 
Pr), with two-dimensional geometrical frustration have 
attracted much attention due to their possible spin-liquid 
state. '^"^ 

Nd3Ga5SiOi4 (NGSO) crystallizes in the trigonal 
space group P321, and Nd'^+ spins {S = 9/2) occupy 
a network of corner-sharing triangles, forming a dis- 
torted kagome-like lattice in the ab plane stacked along 
the c axis, which is topologically equivalent to a per- 
fect one when only the nearest-neighboring interaction 
is considered)^ The magnetic anisotropy associated with 
crystal field effect was confirmed from the magnetic sus- 
ceptibility and nuclear magnetic resonance (NMR) mea- 
surements, which displayed a distinct crossover from the 
easy plane (kagome plane) to the easy axis (c axis) be- 
low 33 Ki^iiS The experimental investigations on the 
ground state of NGSO seemed to be rather confusing. 
The muon spin relaxation (//SR) and neutron scattering 
experiments did not detect any phase transition or long- 
range magnetic order down to milliKelvin temperatures 
in zero fieldi^'^'i^ In some earlier works, the magnetic 
susceptibility showed that the Curie- Weiss temperature 
is about -52 K or -62 K and therefore the frustration 
factor is large (> 1300), indicating a spin-liquid ground 
state with rather strong spin interactions However, 
by measuring the low-lying crystal-field levels of Nd"^"*" 



via the magneto-optical spectroscopy, it was found that 
the exchange interaction of Nd^+ spins is surprisingly 
small, leading to a revised frustration factor of ~ Ijii 
This finding is essentially consistent with the single-ion 
quantum relaxation at low temperatures revealed by the 
neutron scatteringii^ These results ruled out NGSO as a 
spin liquid at least down to 100 mK. So the ground state 
of NGSO remains to be fully understood. In addition, 
the effect of magnetic field on the ground state of NGSO 
has been probed by magnetization and neutron scatter- 
ing measurements At 1.6 K, the magnetization for 
H \\ c and iJ _L c showed a simple spin polarization, but 
the saturated moments are only half of the free Nd'^"'" 
moment, which is still an open question4 d°i^^ 

Low-temperature heat transport has recently been 
found to be an effective way to probe the mag- 
netic excitations and the field-induced quantum phase 
transitions. ^^"^^ In this work, we study the heat trans- 
port properties of NGSO at low temperatures down to 0.3 
K and in magnetic field up to 14 T. The thermal conduc- 
tivity (re) shows a typical phononic heat transport behav- 
ior in zero field and at very low temperatures, while it ex- 
hibits remarkably anisotropic magnetic-field dependence. 
The magnetic field along the c axis causes a low-field dip 
in k{H) isotherms and the high-field k are the same as 
the zero-field values, which can be simply understood as a 
paramagnetic scattering on phonons. Whereas the k{H) 
isotherms in magnetic field along the ab plane show sharp 
and step-like decreases, which indicates a field-induced 
magnetic transition. These transport results point to a 
paramagnetic ground state of NGSO but a field-induced 
partial order for H \\ ab, which is consistent with what 
the low-T specific heat data suggest. 
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II. EXPERIMENTS 

High-quality NGSO single crystals are grown by the 
floating-zone technique in flowing argon and oxygen mix- 
ture with the ratio of 10: To reduce the bubble during 
the growth, the sintered feed rod is firstly pre-melted with 
a rate of 25 mm/h. Then the single-crystal rods (5 mm in 
diameter and 9 mm in length) are successfully obtained 
with a low growth rate of 1 mm/h. It should be noted 
that this fast-scan procedure is very useful to improve 
the quality of NGSO crystals in the aspects of both the 
size and the crystallinityi^i^ The samples for thermal con- 
ductivity measurements are cut into parallelepiped shape 
with a typical size of 3 x 0.6 x 0.15 mm^ along either 
the a axis or the c axis. The thermal conductivity is 
measured at low temperatures down to 0.3 K using a 
conventional steady-state technique, which has been de- 
scribed in details elsewhere j^'i^riS The specific heat is 
measured by the relaxation method in the temperature 
range from 0.4 to 30 K using a commercial physical prop- 
erty measurement system (PPMS, Quantum Design). 



III. RESULTS AND DISCUSSION 

Before presenting the heat transport results of NGSO 
single crystals, we study the low-T specific heat in de- 
tails. Figure 1 shows the temperature dependencies of 
the specific heat in zero field and in magnetic fields for 
H \\ c and _ff _L c. The zero-field specific heat shows a 
upturn at low temperatures below 2 K, which has been 
observed in previous work and is apparently of magnetic 
origin. The distinct peaks develop in various magnetic 
fields both for H \\ c and iJ _L c with the peak posi- 
tion shifting to higher temperatures with increasing field. 
This is a characteristic behavior of the Schottky anomaly, 
which commonly exists in many magnetic materials. The 
phononic specific heat can be estimated by fitting the 
high-T data in zero field using the low-frequency expan- 
sion of the Debye function Cp = ^T^ + /S^T^ + /SyT^^ 
as shown by the dashed lines in Fig. 1. The best fitting 
gives /3 = 1.5 X 10-3 J/K^mol, = 1.29 x 10"*^ J/K^mol, 
and (3j = —1.93 x 10"^ J/K^mol. The magnetic specific 
heat Cm can be obtained by subtracting the lattice con- 
tribution from the total Cp, as shown in Fig. 2. The 
Schottky anomaly in NGSO can be simply interpreted as 
the splitting of a ground-state doublet of Nd'^+ with the 
effective spin S = 1/2. Therefore, the data in magnetic 
fields are tried to fitted by using the two-level Schottky 
anomaly2£iSi 

Csch = Ni]^) (i + eAABT)2' (1) 

where A is the gap value for the Zeeman splitting of the 
Nd'^"'" Kramers doublet. The concentrations of free spins 
is N/R, with R the universal gas constant. 
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FIG. 1: (Color online) Temperature dependencies of the low- 
temperature specific heat of NGSO single crystals in the mag- 
netic field applied along (a) and perpendicular (b) to the c 
axis. The dashed black lines show the fitting of the high-T 
zero-field data by using the Debye model of lattice specific 
heat. 

The best fitting results using Eq. ([l} are also shown in 
Fig. 2. It can be seen that the magnetic specific heat for 
H \\ c can be fitted very well using the two- level Schot- 
tky anomaly, indicating that the spins behave mainly like 
a free system. The fitting parameter A « Q.175H for 

II c, as shown in Fig. 2(c), is in good agreement with 
previous results obtained from inelastic neutron scatter- 
ing spectrumi^ii^ Another fitting parameter N = 22.3 
J/Kmol is nearly 3R, which indicates that almost all the 
Nd'^"'" moments behave like paramagnetic spins in mag- 
netic field along the c axis. For H L c, the similar data 
analysis are also performed and shown in Fig. 2. The 
fitting parameter A ft! 0.125H suggests the anisotropic 
Lande factor. The parameter N — 20 A J/Kmol is rather 
consistent with the value for H \\ c. However, it is no- 
table that there is a pretty larger discrepancy between 
the fitting curves and the specific heat data at very low 
temperatures, particularly for > 1 T, as shown in Fig. 
2(b). This means that at very low temperatures, there 
must be some other kind of magnetic contribution to the 
specific heat. These data suggest that the ground state 
of NGSO is likely paramagnetic rather than spin-liquid 
like. The unexpected phenomenon is that the magnetic 
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FIG. 2: (Color online) Temperature dependencies of mag- 
netic specific heat Cm in field applied along (a) and perpen- 
dicular to (b) the c axis, respectively. The solid lines cor- 
respond to the fitting results by a simple two-level Schottky 
model, (c) The relationship between the fitting parameter 
A and magnetic field. The thin lines show linear dependen- 
cies A = O.I75ff and A = 0.I25H for H \\ c and H ± c, 
respectively. 



field along the ab plane induces some magnetic excita- 
tions or spin fluctuations, which indicate some kind of 
field-induced partial or short-range magnetic order (no 
sign for long-range order in the specific heat). As we can 
see below, such kind of in-plane field induced order has 
remarkable impact on the low-T heat transport. 

Thermal conductivity as a bulk measurement tech- 
nique is a sensitive probe to all elementary excitations 
that carry heat or scatter heat carriers.— Figure 3 shows 
the temperature dependencies of k with heat fiowing 
along the a and c axes in zero field. Both curves show 
phonon peaks at about 7 K, a typical feature of insu- 
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FIG. 3: (Color online) Temperature dependencies of thermal 
conductivity along the a axis and the c axis of NGSO single 
crystals in zero field, respectively. The solid lines are the 
fitting results using the Debye model. The inset shows the 
temperature dependencies of the ratio of the phonon mean 
free path I to the averaged sample width W. 



lating crystals, and exhibit rather weak anisotropy. At 
subKelvin temperatures, a rough T^'^ dependence is ob- 
served, indicating that k{T) behaviors are close to the 
boundary scattering limit. In order to judge whether 
the phonons are free from microscopic scattering at sub- 
Kelvin temperatures, we make an estimation of the mean 
free path I of phonons at low temperaturesJ^ii^ The 
phononic thermal conductivity can be expressed by the 
kinetic formula Kph — ^Cvpl, where C = (3T^ is the 
phonon specific heat at low temperatures and Vp is the 
average velocity. With /3 value obtained from the above 
fitting of the specific-heat data, the average sound ve- 
locity Vp = 2410 m/s and the Debye temperature Qd 
= 310 K are calculated by using the expression &d = 
^(62^)1 and l3 = ^-§§-, where V is the volume of 
crystal, R is the universal gas constant, N is the num- 



ber of molecules per mole, and s is the number atoms 
in each molecule. The inset to Fig. 3 shows the tem- 
perature dependencies of the ratio l/W for Ha and Kc, 
where W is the averaged sample width. It is found that 
both the l/W ratios are rather close to 1 with lowering 
temperature down to 0.3 K, indicating that the phonon 
heat transport of NGSO indeed approaches the boundary 
scattering limit at very low temperatures and the micro- 
scopic scatterings on phonons are nearly wiped out. 

For NGSO, the ^SR and neutron spectroscopy revealed 
that the magnetic fluctuations slow down below 10 K but 
keep persistent down to milliKelvin temperatures in zero 
fleld.^'^^ However, the k.{T) curves do not show any clear 
anomalies around 10 K, which suggests that the magnetic 
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excitations do not couple with phonons. Therefore, it is 
possible to try a more quantitative analysis on the zero- 
field data by using a classical Debye model of phonon 
thermal conductivity^^"— 

in which x = hw/kBT is dimensionless, w is the phonon 
frequency, and t{u}, T) is the mean lifetime of phonon. 
The phonon relaxation is usually defined as 

= Vp/L + Aoj* + BTw3exp(-eD/6r) + t"^^, (3) 

where the four terms represent the phonon scattering 
by the grain boundary, the point defects, the phonon- 
phonon Umklapp scattering, and the resonant scattering, 
respectively. Since the k{T) curves do not show clear 
signature of resonant scattering, the last factor can be 
neglected. The parameters Vp and are obtained from 
the specific-heat data and other parameters L, A, B and 
h are free ones. It is turned out that the zero-field ther- 
mal conductivity below 25 K is well simulated as the blue 
lines shown in Fig. 3. The best fitting parameters are 
L = 3.7 X 10-4 m, A = 8.15 x IQ-*! s^, B = 9.8 x 10-^9 
K-is2, b = 8.7, and L = 4.8 x 10"^ va, A = 5.5 x IQ-^i 
s3, B = 5.0 X 10-29 K-is2, b = 9.0 for and Kc, re- 
spectively. The parameters L are found to be in good 
consistent with the averaged sample widths of these two 
samples, which are 3.70 x 10"^ m and 3.95 x 10"^ 
m (kc), respectively. Note that the failure of fitting at 
T > 25 K can be related to either the magnetic fluctua- 
tions or the contributions of optical phononsi^ However, 
these complexities at relatively high temperatures are not 
involved in the following discussion about the magnetic- 
field effect on low-T heat transport. 

Detailed magnetic-field dependencies of the low-T 
thermal conductivity can provide useful information on 
the nature of the ground state. Figure 4(a) shows the 
detailed magnetic field dependencies of Ka for H \\ c. 
The n{H) are strongly suppressed down to about 50% 
at low magnetic field and are recovered at high mag- 
netic field, yielding a dip-like feature. Note that the 
position of the minimum moves to high magnetic field 
as increasing temperature. Based on the Schottky con- 
tribution to the low-T specific heat, the dip- like behav- 
iors of k{H) can be attributed to the phonon scatter- 
ing by the paramagnetic moments, that is, the field- 
induced splitting of Nd'^"'" Kramers doublet can cause 
a resonant scattering on phonons, as we previously ob- 
served in Pri.3Lao.7Cu04 and GdBaCoaOs+x'^'^e The 
Zeeman splitting of the spin doublet is known from the 
above fitting of Schottky anomaly for H \\ c, which gives 
A(iJ) « 0.175iJ. We can simulate the magnetic-field de- 
pendencies of k{H) by using the Debye model for phonon 
thermal conductivity,— as shown in Fig. 4(b). One can 
see that most of k{H) isotherms are in good agreement 
with the calculated ones except for some deviations for 
the high-T data, which might be due to the excitations 
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FIG. 4: (Color online) (a) Magnetic-field dependencies of 
thermal conductivity of NGSO single crystal at low tempera- 
tures with magnetic field along the c axis, (b) The calculation 
results of k{H) isotherms using the Debye model. 



of the higher-energy crystal field levels. Therefore, the 
c-axis magnetic field seems to affect the heat transport 
mainly by introducing some paramagnetic scattering on 
phonons. This is in good correspondence with what the 
specific-heat data suggest. It is worthy of pointing out 
that the Kc{H) isotherms for || c are the same as those 
of Ka{H) shown in Fig. 4. 

The above paramagnetic scattering effect is known to 
be qualitatively isotropic for magnetic field along differ- 
ent directions 1^1^ However, the heat transport behaves 
very differently if the magnetic field is applied in the ah 
plane. Figure 5 shows the magnetic-field dependencies 
of both Ka and Kc for H J- c. It is easily found that 
the k{H) behaviors are essentially the same for different 
directions of heat current. All the k(H) isotherms show 
significant suppression down to 10-20% of the zero-field 
value. With lowering temperature, the decrease of k be- 
comes sharper at low fields. If a dip- like k{H) term from 
the paramagnetic scattering effect is also included, the 
field dependence of Fig. 5 apparently has another con- 
tribution. This can be clarified by a simple analysis. As 
some examples shown in Fig. 6 (similar analysis can be 
done for other data), the changes of Kq induced by para- 
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FIG. 5: (Color online) Magnetic-field dependencies of thermal 
conductivities K,a (a) and Kc (b) of NGSO single crystals at 
low temperatures with magnetic field perpendicular to the c 
axis. 



magnetic scattering are calculated in the same way as 
that for H \\ c, considering A(iJ) w 0.125iJ. Then sub- 
tracting this part of change from the raw data of Ka{H), 
the other contribution to the field dependence can be 
roughly estimated. It can be seen that at low tempera- 
tures, this additional term is actually a step-like suppres- 
sion, which suggests that the a&-plane field induces some 
magnetic transition. This is compatible with what the 
specific-heat data signify, as discussed above. Although 
the quantitative description of the phonon scattering by 
the magnetic excitations in such a broad temperature re- 
gion is not available now, the qualitative understanding 
can be easily figured out. Apparently, the field-induced 
partial or short-range order causes strong dynamic spin 
fluctuations that can strongly scatter phonons. 

It should be pointed out that the in-plane field induced 
partial order, suggested by the specific heat and thermal 
conductivity data in the present work, however, has not 
been evidenced in the earlier works. It is known that 
the low-T magnetization showed a simple spin polariza- 
tion above - 2 T for both H \\ c and H _L and 
the neutron scattering also confirmed this phenomenon 
for 7? II But the moments of the high-field M{H) 
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FIG. 6: (Color online) Magnetic-field-induced changes of Ka 
at 0.36 K (a) and 0.5 K (b) for _L c. A«:(i7) is defined 
as K,{H) — k(0). The solid symbols indicate the experimental 
data. The solid blue lines represent the calculated changes 
of K, by the paramagnetic scattering, denoted as AnpAf. The 
open symbols show the results of subtracting the paramag- 
netic scattering term from the raw data, denoted as Ahfio, 
which is discussed to be related to some kind of field-induced 
magnetic order. 



plateau are only half of the free Nd'^"'" moment, which is 
an open question for the existing theories on the kagome 
antifcrromagneti^iiSii^ There have been no other experi- 
mental investigation about the details for H ± c. Here, 
we propose one possibility that when the Nd^^^ moments 
are gradually aligned along the in-plane field, the c-axis 
component of Nd"^+ moments may form some kind of AF 
ordering. They could be anti-parallel either for the neigh- 
boring sites or for the neighboring layers. This picture 
is not inconsistent with the experimental results that the 
high- field M{H) plateau is lower for H J- c than that 
for H II c— Nevertheless, the validity of this canted 
AF order and the underlying mechanism, which should 
be related to the spin anisotropy, the crystal field effect, 
and the spin frustration, calls for a deeper investigation. 
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IV. CONCLUSION 

In summary, the temperature and magnetic field de- 
pendencies of the thermal conductivity and specific heat 
are studied for NGSO single crystals. The zero-field ther- 
mal conductivity shows a purely phononic heat trans- 
port behavior at low temperatures, without showing the 
strong coupling between phonons and magnetic excita- 
tions. A dip-like field dependence of k for magnetic field 
along the c axis can be explained by a simple param- 
agnetic scattering on phonons; whereas a step-like field 
dependence of k{H) for iJ ± c indicates a field-induced 
magnetic transition at low temperatures. These behav- 
iors are consistent with the field dependencies of low-T 
specific heat. The present results demonstrate a param- 



agnetic ground state rather than the spin liquid. How- 
ever, the partial magnetic order induced by the field along 
the ab plane seems to be quite unusual and remains to 
be further investigated both experimentally and theoret- 
ically. 
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